简化 FMM 算法 在 束 流 剖面 测量 中 的 应 用 
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用 荧光 靶 系 统 测量 束 流 横向 剖面 ， 为 进一步 提高 系统 测量 的 准确 度 ， 本 文 提出 
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析 了 算法 处 理 前 后 的 模拟 图 像 束 流 剖 面 参数 ， 仿 真 结 果 表 明 经 过 处 理 后 的 模拟 图 像 ， 其 结构 相似 
Index measure, SSIM) 更 接近 1、 均 方 误差 (Mean Square Error, MSE) 


医用 重 离子 加 速 器 HIMM (Heavy Ion Medical Machine) 现场 采集 的 荧光 靶 图 像 
进行 验证 。 结 果 表 明 ，HIMM 现场 采集 的 图 像 经 处 理 后 ， 剖 面 测量 结果 接近 评估 值 。 本 研究 结果 


中 图 分 类 号 : TL56, TP751 文献 标志 码 : A 


DOI: 


HIMM 是 由 中 国 科学 院 近代 物理 研究 所 和 兰州 科 近 泰 基 新 技术 有 限 责 任 公司 研发 制 
的 大 型 放疗 设备 03。 在 调 束 过 程 中 ， 调 束 人 员 为 及 时 了 解 束 流 状态 信息 ， 多 采用 简 
系统 获取 束 流 横向 剖面 参数 。 荧 光 靶 靶 片 在 安装 前 会 采用 石墨 烧 制 技术 刻 制 网 格 状 


1 引言 
医用 重 离子 加 速 器 

造 、 拥 有 自主 知识 产权 

单 、 直 接 的 荧光 靶 测 量 

EAC, 

通过 网 格 标定 线 获取 束 


满足 实时 监测 束 流 拌 动 情况 的 需求 ， 同 时 标定 线 会 给 束 流 训 面 参数 计算 引入 误差 S14， 因此 ， 为 实现 准确 、 


网 格 尺寸 为 5mm*5mm， 束 流 打 在 训 面 上 发 出 可 见 光 ， 并 由 相机 采集 传输 至 处 理 器 ， 调 束 人 员 
流 位 置 和 剖面 信息 。 相 机 帧 率 为 25 帧 每 秒 ， 然 而 ， 上 位 机 程序 处 理 时 间 较 长 ， 无 法 


快速 的 参数 计算 ， 同 时 考虑 FPGA 算法 实现 难 易 程度 ， 本 文采 用 简化 FMM 算法 填充 标定 线 并 设计 了 计算 光 
斑 中 心 位 置 和 标准 差 参 数 的 图 像 处 理 流 程 ， 实 现 了 标定 线 填充 和 束 流 齐 面 的 测量 ， 可 为 调 束 人 员 提 高 调 束 
测量 效率 。 
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为 提高 系统 测量 的 准确 度 ， 本 文 在 该 系统 上 设计 图 像 处 理 算法 中， 并 验证 算法 的 可 用 性 。 
3 ”标定 线 图 像 填 充 

HRA BA CE VARIN Se A LL As BY EERE BP RA J, MELE Pa TD Ge, ee 
的 存在 导致 计算 给 出 的 剖面 参数 与 实际 参数 存在 一 定 误 差 ， 因 此 ， 需 要 利用 数字 图 像 修 复 (digital image 
inpainting ) 技术 填充 标定 线 ， 提 高 测量 结果 的 准确 性 。 数 字 图 像 修 复 是 图 像 处 理 领域 的 重要 分 支 之 一 ， 这 
一 概念 由 Bertalmio 等 于 2000 年 明确 提出 501， 是 指 利 用 已 知 信息 对 图 像 中 缺损 区 域 进 行 填充 的 一 种 技术 。 
在 传统 小 尺度 缺损 修复 方法 中 ，Telea 提出 的 FMM 相对 其 他 算法 实现 相对 简单 、 修 复 速度 较 快 0223]， 东 
南大 学 的 闵 溪 青 在 此 基础 上 ， 提 出 了 一 种 简化 的 FMM 算法 0， 与 FMM 相 比 ， 其 保持 了 修复 效果 的 同时 
减少 了 运算 量 ， 且 算法 实现 更 为 简单 。 
3.1 简化 FMM 算法 基本 原理 
简化 FMM 算法 的 原理 图 如 图 2 所 示 ， 其 中 Unknown 为 待 修复 区 域内 部 ， 该 部 分 像素 点 的 像素 值 未 

知 ; Known 为 待 修复 区 域外 部 ， 该 部 分 像素 点 的 像素 值 已 知 ，q 点 为 待 修复 点 ， 标 记 为 Unknown; N(q) 
AVA q 点 为 中 心 的 4 邻 域 区 域 。 简 化 FMM 算法 的 修复 流程 如 下 : 


Known wx, 


UnkKonwn 


图 2 简化 FMM 算法 原理 图 

1) 遍历 掩 码 图 像 〈 像 素 值 只 有 0 和 255) ， 像 素 值 为 0 的 点 待 处 理 图 像 相 应 位 置 标记 为 Known 的 
像素 ， 像 素 值 为 255 的 点 则 标记 为 Unknown 的 像素 。 

2) 遍历 待 修复 区 域 ， 对 于 Unknown 上 的 任 一 点 q， 若 点 q 的 四 邻 域 像素 点 均 为 Unknown 标记 ， 则 
跳 过 。 若 p 的 四 邻 域 像 素 值 存 在 n 个 标记 为 Known 的 点 ， 根 据 式 1 填充 点 q， 填 充 后 将 q 点 标记 为 
Known. 

3) 更 新 待 修复 区 域 并 返回 步骤 2) ， 直 至 待 修复 区 域 为 空 。 


> 1 


I(q)= ieN(q) (1) 
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3.2 修复 前 图 像 预 处 理 
图 像 预 处 理 的 主要 目的 是 提高 输入 图 像 质量 和 分 割 标定 线 生 成 掩 码 图 。 本 文 的 预 处 理 方法 主要 包括 
灰 度 化 、 滤 波 降 品 、 透 视 变换 、 自 适应 闵 值 分 割 、 腐 蚀 膨 胀 等。 图像 处 理 流程 图 如 图 3 所 示 。 
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变 的 荧光 靶 图 像 和 标定 线 掩 码 图 像 ， 其 作为 简化 FMM BARERA RR, h RAE 2a RGE 


图 像 。 
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预 处 理 流程 


， 灰 度 化 操作 可 以 减少 数据 量 ， 中 值 滤波 可 以 去 除 椒盐 噪声 的 同时 保护 图 像 边 缘 信 ， 
于 系统 真空 管道 的 机 械 结 构 ， 荧 光 刘 探 针 在 安装 时 不 
现场 安装 时 的 误差 ， 从 而 相机 采集 到 的 图 像 发 
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首先 计算 出 3x3 


3.3 ”模拟 荧光 靶 图 像 填充 效果 


xl 


采 月 


光 靶 图像 ， 并 将 该 图 


H Matlab 语言 编写 代码 4 


成 640*480 大 小 的 光斑 图 
中 心 坐标 为 250; 水 平方 向 标准 差 为 30， 中 心 坐标 为 300。 在 光斑 图 像 


的 标定 线 ， 生 成 掩 码 
可 使 掩 码 区 域 增 大 。 预 处 理 结 
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+ 滤波 降 噪 一 。 透视 变换 
”简化 FMM _ 
修复 
参数 计算 
图 3 图 像 处 理 流程 图 


Se 
能 将 相机 、 


个 变换 矩阵 ， 利 用 变换 矩阵 将 一 个 像 平 
邻 域 块 的 算术 平均 值 ， 再 乘 以 一 个 缩小 
图 像 ， 腐蚀 操作 可 以 消除 掩 码 图 
J. 
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束 得 到 消除 噪声 日 


Bam 


像 进行 算法 测试 ， 光 斑 的 垂直 方向 标准 差 为 25， 


加 入 噪声 和 标定 线 生 成 模拟 的 区 


像 作为 透视 变换 后 的 荧光 靶 图 像 ， 模 拟 图 像 如 图 4 所 示 。 
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通过 简化 FMM 算法 修复 后 的 结果 如 图 


5 所 示 。 
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填充 模拟 荧光 靶 图 像 效 果 


引入 结构 相似 性 和 均 方 误差 客观 指标 衡量 填充 效果 053。 结 构 性 相似 SSIM， 是 一 种 度量 两 幅 图 像 相似 
度 的 指标 。SSIM 算法 主要 比较 两 图 像 的 亮度 、 对 比 度 、 结 构 三 要 素 ， 然 后 对 三 要 素 加 权 乘 积 表示 。 式 2 
为 SSM 的 简化 公式 : 


(2) 
SSIM (x, y) = (24.4, taea, ra) 
(u +u, +eco +o, +c) 
1 N 
Oy =F TÀ (x; = (V: = 4, ) 
CHAMIAR x My, AP, HARRNBA PH, o 为 图 像 的 像素 方差 ，N 为 图 像 像 素 总 数 ， 


关 和 用 分 别 代表 图 像 该 点 的 像素 值 ，c 和 名 是 保证 公式 稳定 性 的 常数 ，c = (KiL)? c,=(K,LY, LX 


图 像 的 灰 度 级 数 ， 对 于 8bit 的 灰 度 图 , 工 取 255$， 及 , 常 取 0.01， 天 ， 常 取 0.03。SSIM 的 范围 是 0 到 1， 两 
图 像 越 相似 其 SSIM 值 越 接 近 1， 当 SSIM 值 为 1 时 ， 表 示 两 图 像 相同 。 
均 方 误差 MSE， 是 用 来 计算 两 幅 图 像 每 个 点 位 置 上 的 像素 差异 。 其 公式 如 式 3: 


m-l n-l 
MSE = 一 》 FG, -KG AT G) 
i=0 j=0 
已 知 有 两 幅 图 像 T 和 KK， 其 中 ，m、n 为 对 应 的 像素 总 数 ，7(i, j) 和 天 (记分 别 代表 两 幅 图 在 (7 
处 的 像素 值 。 在 判断 填充 质量 时 ，MSE 值 越 小 代表 修复 效果 越 好 。 表 1 为 采用 简化 FMM 算法 前 后 的 指 
标 对 比 ， 由 表 1 可 知 ， 填 充 后 的 SSIM 值 更 接近 1、MSE 更 小 ， 表 示 填 充 后 的 荧光 车 图 像 与 原 光斑 图 像 更 
加 相似 。 


表 1 简化 FMM 算法 前 后 指标 对 比 


SSIM MSE 
去 标定 线 前 0.946 114.872 
去 标定 线 后 0.996 0.858 


图 6、 图 7 分 别 为 标定 线 填充 前 后 水 平方 向 和 垂直 方向 的 分 布 ， 结 合 图 5-7 观察 明显 得 出 ， 填 充 后 的 
奕 光 靶 图 像 与 原 光斑 图 像 分 布 更 近似 。 
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图 6 填充 前 后 水 平方 向 分 布 折线 图 
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图 7 填充 前 后 垂直 方向 分 布 折线 
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4 ， 束 流 剖面 参数 分 析 
束 流 的 中 心 位 置 和 标准 差 是 束 流 剖 面 测量 的 重要 参数 ncia。 对 于 束 流 而 言 ， 其 理想 的 光斑 图 像 的 光 
强 分 布 应 符合 二 维 高 斯 分 布 分 布 函数 为 


He) =e EM +07) | (4) 
o 


x 了 


Er 


1: 了 (Xx,y) 为 束 流 截 面 (x,y) 处 的 光 强 ， 对 应 图 像 坐标 (x, y) 处 的 像素 值 ， 4 为 束 流 截 面 的 光 强 幅 值 ; 
(Xo Vo) 为 光斑 中 心 位 置 ， G0 ，a， 分 别 为 水 平方 向 和 垂直 方向 的 标准 差 ; C 为 偏 置 常数 。 依 据 二 维 高 斯 
分 布 的 边缘 概率 密度 特性 ， 对 束 流 图 像 的 水 平和 垂直 方向 分 别 投影 获得 两 方向 的 一 维 高 斯 分 布 。 例 如 水 
平方 向 的 分 布 函数 g(xX) 为 


g(x) => (x,y) (5) 


同时 式 5 满足 一 维 高 斯 分 布 ， 分 布 函数 为 


由 
g(x) = A, -exp)— +C, (6) 


式 中 : 4 为 水 平方 向 的 幅 值 ; C, 为 水 平方 向 的 偏 置 。LM 〈 列 文 伯 格 - 马 夸 尔 特 ，Levenberg- Marquardt) 
算法 L920 是 一 种 非 线 性 最 小 二 乘 算法 ， 常 用 的 软件 ， 如 : Origin, Matlab. LabVIEW 等 都 采用 LM 算法 进 
行 非 线性 拟 合 ， 且 LM 算法 提出 一 种 信赖 域 思想 ， 控 制 了 迭代 步 长 ， 缩 短 了 迭代 时 间 。 本 文选 择 LM 算 
去 进行 高 斯 拟 合 求 参 。 在 实际 测试 过 程 中 ， 由 于 LM 算法 是 局 部 最 优 算法 ， 设 置 的 初 值 对 拟 合 结果 影响 
很 大 ， 常 常 初 值 设置 不 合适 ， 导 致 无 法 拟 合 出 正确 结果 。 

束 流 的 半 高 全 宽 也 是 束 流 剖 面 的 重要 参数 ， 当 束 流 的 形状 不 规则 时 ， 通 常 实 验 人 员 以 束 流 的 半 高 全 
宽 作为 束 流 的 剖面 参数 。 经 过 实验 测试 ， 以 束 流 的 半 高 全 宽 作为 拟 合 的 初 值 可 以 得 到 正确 的 拟 合 结 
在 已 知 束 流 某 方向 的 分 布 函 数 g(X) 后 ， 可 以 在 该 函数 内 找到 两 点 x 和 x,， 使 式 7 成 立 。 


ee 


= 


g(x) ¥ g(x) ae Zrin (7) 


其 中 Zna 和 gi 分 别 是 分 布 函数 g(x) 的 最 大 值 和 最 小 值 ， 束 流 在 该 方向 的 半 高 全 宽 OY 


6, =|, -| (8) 


同 理 可 以 得 到 垂直 方向 的 束 流 半 高 全 宽 。 束 流 该 方向 的 标准 差 o 与 0. 的 关系 为 


5, = 2.3550 (9) 


g 的 坐标 索引 为 束 流 的 中 心 位 置 ，o 为 束 流 的 标准 差 。 表 2， 表 3 分 别 给 出 了 模拟 荧光 靶 图 像 使 
用 简化 FMM 算法 填充 标定 线 前 后 剖面 参数 的 计算 结果 和 误差 数值 比较 ， 其 中 中 心 位 置 取 四 使 五 入 后 的 结 
果 ， 结 果 显 示 填 充 后 的 计算 结果 更 接近 原 光 斑 分 布 ， 说 明 本 方法 获得 的 参数 计算 结果 具有 较 高 可 信 度 。 


表 2 简化 FMM 算法 前 后 参数 对 比 


算法 中 心 位 置 ( 单 ”垂直 方向 标准 ”水 平方 向 标准 
位 : 像素 ) 差 ( 单 位: 像 ” 差 (单位 : 像 
素 ) 素 ) 
模拟 值 (250,300) 25 30 
填充 前 (250,300) 26.484 32.242 
填充 后 (250,300) 24.764 29.711 


表 3 简化 FMM 算法 前 后 误差 对 比 


算法 MARA ( 单 ”垂直 方向 标准 ”水 平方 向 标准 
位 : %) ZRA (Ff 差 误差 ( 单 
位 : %) 位 : %) 
填充 前 0 5.9 75 
填充 后 0 0.9 0.9 


5 HIMM 荧光 靶 图 像 效果 验证 


HIMM 现场 采集 到 的 荧光 靶 光 斑 图 像 (如 图 8) 分 辨 率 为 1920x1080， 由 于 安装 角度 原因 ， 靶 片 网 格 
由 正方 形变 成 了 长 方形 ， 选 取 束 流 中 心 4X6 个 网 格 区 域 的 四 顶点 进行 透视 变换 ， 透 视 变 换 矩 阵 为 


[0.5134, -0.03305, -220.8; -0.01668, 0.8674, -435.53; 0.0000142，-0.000194，1]， 变 换 后 图 像 分 辨 率 为 
350x234， 如 图 9 所 示 。 
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图 8 HIMM 现场 采集 的 效 光 靶 光 斑 图 


图 9 3B LAR a EAE EE 
选取 中 心 2x4 网 格 区 域 交点 间 连 线 ， 以 水 平方 向 为 例 ， 计 算 标定 线 间 宽度 为 59.3 土 0.5 个 像素 ， 由 于 


标定 线 间 宽度 计算 存在 误差 ， 导 致 计 算 标定 值 时 的 误差 传递 ， 误 差 传递 公式 如 式 10: 


10 
2 2 二 2, (omy a 


式 中 : mp 为 逐 像素 标定 值 (单位 mm/pixel), WA 5 为 标定 线 间 宽度 和 高 度 差 (单位 mm), Op 为 逐 


像素 标定 值 的 标准 差 〈 单 位 mm/pixel) ， 克 为 两 标定 线 间 水 平 宽度 的 均值 。 计 算得 水 平方 向 逐 像素 的 标 
定 值 为 0.084 土 0.000318 mm/pixel。 同 理 ， 计 算得 两 标定 线 间 高 度 为 60.6 土 0.7 个 像素 ， 逐 像素 标定 值 为 
0.083 +0.000426 mmyApixel。 透 视 变 换 后 束 流 中 心 位 置 坐标 (162, 115) ; 水 平方 向 标准 差 为 118.36 个 像 
素 即 9.94 士 0.04mm， 垂 直方 向 标准 差 为 30.66 个 像素 即 2.54 士 0.01mm 。 透 视 变换 后 不 使 用 简化 FMM 算 
法 填充 标定 线 ， 计 算 的 光斑 参数 为 ， 中 心 位 置 坐标 (161.86, 115.96) ; 束 流水 平方 向 标准 差 为 99.854 个 
GAM 8.32mm， 垂 直方 向 标准 差 为 29.02 个 像素 即 2.42mm。 透 视 变 换 后 采用 自 适 应 阔 值 分 割 获取 标 定 线 
掩 码 图 像 ， 如 图 10 左 图 所 示 。 简 化 FMM 算法 处 理 后 的 结果 如 图 10 右 图 所 示 。 填 充 标 定 线 后 计算 的 光斑 
参数 为 :中心 位 置 坐 标 (161.92, 115.09) ; 束 流水 平方 向 标准 差 为 114.89 个 像素 即 958mm, HAT E 
标准 差 为 30.23 个 像素 即 2.52mm。 
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图 10 掩 码 图 E) 、 fide EMM sb 后 结果 图 A) 


综 上 ， 简 化 FMM 算法 填充 标定 线 前 后 中 心 位 置 计算 结果 近似 ， 对 比 结果 表示 ， 填 充 后 标准 差 的 计算 
结果 更 加 准确 ; 填充 前 由 于 标定 线 的 影响 ， 标 准 差 水 平方 向 的 计算 误差 约 为 19.47%， 垂 直方 向 的 计算 误 
差 约 为 4.96%; 填充 后 标准 差 水 平方 向 的 计算 误差 约 为 3.76%， 垂 直方 向 的 计算 误差 约 为 0.79%， 计 算 误 
差 约 减少 到 填充 前 的 14， 说 明 简 化 FMM 算法 可 以 有 效 提 高 剖面 测量 精度 。 但 由 图 11 可 见 ， 简 化 FMM 
算法 填充 标定 线 的 结果 并 不 符合 人 眼 主观 要 求 ， 这 是 由 于 本 文 测试 采用 的 HIMM 现场 荧光 靶 图 像 是 在 开 
J 条 件 下 拍摄 的 ， 荧 光 靶 靶 面 背景 灰 度 值 在 110-180 之 间 ， 因 此 自 适应 闵 值 分 割 算法 会 将 非 光斑 周围 的 标 


定 线 提取 出 来 ， 少 部 分 没有 提取 出 来 的 标定 线 周 围 像 素 点 会 在 计算 填充 值 时 与 背景 像素 进行 加 和 平均 ， 
从 而 影响 填充 结果 。 
6 结论 

本 文 根 据 项 目 实 际 需求 ， 将 简化 FMM 算法 应 用 在 菊 光 靶 束 流 剖 面 测量 中 。 模 拟 和 现场 测试 结果 表明 ， 
使 用 简化 FMM 算法 去 除 标定 线 后 的 欧 光 刘 图 像 计算 结果 与 实际 结果 更 接近 ， 计 算 误 差 减 少 了 近 34， 且 
简化 FMM 算法 易于 实现 ， 可 以 提高 系统 测量 的 准确 度 。 
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Application of simplified FMM algorithm in beam profile 


measurement 
ZHANG Hongbo!?, WEI Xiaojuan!*, LI Min?3*, MAO Ruishi?*, TANG Kai*, FENG Yongchun? 
(1.Northwest Minzu University, Lanzhou 730030, China ; 
2.Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China 
3. University of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: The fluorescence screen system is used to measure beam transverse profile in the accelerator. In order to further improve 
the measurement accuracy of the system, this paper proposes a method, filling the orientation lines of fluorescence screen with a 
simplified Fast Matching Method (FMM), so that the filled fluorescence screen images can reflect the beam profile information more 
truly. Firstly, the study uses Matlab to generate the simulated fluorescence screen images, and comparative analysis the beam profile 
parameters of the analog image before and after algorithm processing. The simulation results show that the Structural Similarity 
Index Measure (SSIM) of the processed simulated image is closer to 1, and the Mean Square Error (MSE) is smaller. Finally, the 
fluorescence screen images collected by the medical Heavy Ion Medical Machine (HIMM) were used for verification. After 
processing the fluorescence screen images collected by HIMM, the profile measurement results are close to the evaluation value. The 
results of this study provide a basis for the FPGA algorithm acceleration of the subsequent beam transverse profile parameter 
algorithm. 

Key words: beam profile measurement; simplified FMM algorithm; fluorescence screen; orientation line filling; image processing 
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